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Excitations from a strongly frustrated system, the kagome ice state of the spin ice Dy2Ti2 07 
under magnetic fields along a [111] direction, have been studied. They are theoretically proposed 
to be regarded as magnetic monopoles. Neutron scattering measurements of spin correlations 
show that close to the critical point the monopoles are fluctuating between high- and low- 
density states, supporting that the magnetic Coulomb force acts between them. Specific heat 
measurements show that monopole-pair creation obeys an Arrhenius law, indicating that the 
density of monopoles can be controlled by temperature and magnetic field. 
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Since the quantum mechanical hypothesis of the ex- 
istence of magnetic monopoles proposed by Dirac/'^^ 
many experimental searches have been performed, rang- 
ing from a monopole search in rocks of the moon 
to experiments using high energy accelerators.^^ But 
none of them was successful, and the monopole is an 
open question in high energy physics. Recently, theo- 
retical attention has turned to condensed matter sys- 
tems where tractable analogs of magnetic monopoles 
might be found, and one prediction^^ is for an emer- 
gent elementary excitation in the spin ice^' compound 
DysTisOr. 

In solid water, the protons are disordered even at ab- 
solute zero temperature and thus retain finite entropy,^^ 
and spin ice exhibits the same type of disordered ground 
states.^' The Dy spins occupy a cubic pyrochlore lat- 
tice, which is a corner sharing network of tetrahedra 
(Fig. 1(a)). Each spin is parallel to a local [111] easy 
axis, and interacts with neighboring spins via an effec- 
tive ferromagnetic coupling. This brings about a geomet- 
rical frustration where the lowest energy spin configura- 
tions on each tetrahedron follow the ice rule, "2-in, 2-out" 
structure, and the ground states of the entire tetrahedral 
network are macroscopically degenerate in the same way 
as the disordered protons in water ice.^'^^^ In addition 
to this remarkable observation, there is the more intrigu- 
ing possibility^^ that the excitations from these highly 
degenerate ground states are topological in nature and 
mathematically equivalent to magnetic monopoles. 

The macroscopic degeneracy of the spin ice state can 
be partly lifted by applying a small magnetic field along 
a [111] direction. Along this direction the pyrochlore 



lattice consists of a stacking of triangular and kagome 
lattices (Fig. 1(a)). In this field-induced ground state, the 
spins on the triangular lattices are parallel to the field 
and consequently drop out of the problem, while those 
on the kagome lattices retain disorder under the same ice 
rules, only with a smaller zero-point entropy. This is 
referred to as the kagome ice state-^-^'-^^ -^^^ (Fig. 1(b)). 

In Fig. 1(c) we illustrate creation of a magnetic 
monopole and antimonopole pair in the kagome ice state. 
An excitation is generated by flipping a spin on the 
kagome lattice, which results in ice-rule-breaking "3- 
in, 1-out" (magnetic monopole) and "1-in, 3-out" (anti- 
monopole) tetrahedral neighbors. From the viewpoint 
of the dumbbell model, where a magnetic moment is 
replaced by a pair of magnetic charges, the ice-rule- 
breaking tetrahedra simulate magnetic monopoles, with 
net positive and negative charges sitting on the centers of 
tetrahedra. The monopoles should interact via the mag- 
netic Coulomb force,^^ which is brought about by the 
dipolar interaction^^^ between spins in Dy2Ti207. They 
can move and separate by consecutively flipping spins, 
but are confined to the two-dimensional kagome layer 
(e.g. Fig. 1(d)). This possibility of separating the local 
excitation into its constituent parts is a novel fraction- 
alization in a frustrated system in two or three dimen- 
sions,^' -'^^^ and enables many new aspects of these emer- 
gent excitations to be studied experimentally, such as 
pair creation and interaction, individual motion, currents 
of monopoles, correlations and cooperative phenomena. 
In the present study, inspired by the theoretical predic- 
tion of the monopoles, we have investigated two aspects 
of magnetic monopoles in spin ice using direct neutron 
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Fig. 1. (Color) Magnetic moments of Dy2Ti2 07 reside on the py- 
rochlore lattice. At low temperatures, four magnetic moments 
on each tetrahedron obey the ice rule (2-in, 2-out). The resulting 
spin ice state is shown in (a). The pyrochlore lattice consists of 
stacked triangular and kagome lattices, shown by green and blue 
lines, respectively, along a [111] direction, (b) Under small [111] 
magnetic fields, spins on the kagome lattice remain in the dis- 
ordered kagome ice state. ■'^■'^^ (c) An excited state is induced by 
flipping a spin from (b), enclosed by a dashed circle, where neigh- 
boring tetrahedra have 3-in, 1-out and 1-in, 3-out configurations. 
These ice-rule-breaking tetrahedra are represented by magnetic 
monopoles with opposite charges depicted by spheres, (d) By 
consecutively flipping two spins from (c), the monopoles are frac- 
tionalized. (e) As the magnetic fleld is increased, H ^ He, spins 
realize a fully ordered, staggered arrangement of monopoles. 



scattering techniques and thermodynamic specific heat 
measurements. 

Single crystals of Dy2Ti207 were prepared by the 
fioating-zone method. Specific heat was measured by a 
quasi-adiabatic method. Neutron scattering experiments 
on a single crystal under a [111] field were performed on 
the triple-axis spectrometers BT-9 at the NIST Center 
for Neutron Research and the ISSP-GPTAS at the Japan 
Atomic Energy Agency. The sample was mounted in di- 
lution refrigerators so as to measure the scattering plane 
perpendicular to the [111] direction. 

A straightforward signature of monopole-pair cre- 
ation is an Arrhenius law in the temperature (T) depen- 
dence of the specific heat (C), C{T) oc exp(— A£^//cb^), 
where AE is a field {H) dependent creation energy. One 
can simply expect = Eq — jaH owing to the Zeeman 
effect. Figure 2 shows the measured C{T) of Dy2Ti207 
under a [111] applied field as a function of 1/T. The 
Arrhenius law is clearly seen at low temperatures, indi- 
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Fig. 2. (Color) Speciflc heat under [111] flelds is plotted as a func- 
tion of 1/T. In the intermediate temperature range these data 
are well represented by the Arrhenius law denoted by solid lines. 
The inset shows the fleld dependence of the activation energy. 



eating that monopole-antimonopole pairs are thermally 
activated from the ground state. We remark that all the 
measurements were performed under field cooling condi- 
tions, which are important to avoid complications due to 
spin freezing^' among the ground state manifolds, 
whose degeneracy are slightly lifted. We think that the 
deviation from the Arrhenius law at the lowest temper- 
atures is attributable to these perturbative effects. The 
observed activation energy AE depends linearly on H 
(Fig. 2 inset) between 0.2 and 0.9 T, i.e., in the kagome 
ice state. The deviation from linearity below 0.2 T (spin 
ice regime) suggests that the nearest-neighbor effective 
bond energy Jeff slightly changes between the two states. 
The zero-field value AE{H = 0) = 3.5 K reasonably 
agrees with an estimation AE{H = 0) = 4Jeff = 4.5 
K using Jeff in ref. 16. The observed Arrhenius law of 
C(T)^ which is attributable to variation of density of the 
monopole pairs, implies that the number of monopoles 
can be tuned by changing T and H. 

A microscopic experimental method of observing 
monopoles is magnetic neutron scattering. One challenge 
to the experiments is to distinguish the relatively weak 
scattering from the small number of monopoles from the 
very strong magnetic scattering^' of the ground states. 
A theoretical idea^^ which is helpful for identifying the 
monopole scattering is that the [111] field acts as chem- 
ical potential of the monopoles, enabling us to control 
their density as shown by the present specific heat mea- 
surements. As the field is increased, the kagome ice state 
with low monopole density changes continuously to the 
maximum density state, the staggered monopole state 
(Fig. 1(e)), where all spin configurations become "3-in, 
1-out" or "1-in, 3-out" to minimize the Zeeman energy. 

For the present neutron scattering experiments, the 
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Fig. 3. (Color) The magnetic Bragg intensity at T = Tc +0.05 K 
is plotted as a function of the [111] field. The open squares and 
dashed curves represent the measurements at (220) and corre- 
sponding MC simulations, respectively. The dot-dashed curve is 
the density of positively charged monopoles obtained by the sim- 
ulation. The inset shows the HT phase diagram under the [111] 
field. The solid line represents the first-order phase transition 
with the critical point shown by an open circle. -"^^^ The dashed 
lines are crossovers. -"^^^ The intensity maps shown in Fig. 4 were 
measured at the two points depicted by open squares. 



best temperature and field region to observe monopoles 
in Dy2Ti207 is close to the liquid-gas type critical 
point^^^ (Tc,i^c) (Fig. 3 inset). In the monopole pic- 
ture, where they are interacting via the magnetic 
Coulomb force, the first-order phase transition^^^ is as- 
cribed to phase separation between high- and low-density 
states. We naturally anticipate that neutron scattering 
close to the critical point is a superposition of the scatter- 
ing pattern by the (low-density) kagome ice state^^^ and 
that by high-density monopoles, which is diffuse scatter- 
ing around magnetic Bragg reflections, i.e., ferromagnetic 
fluctuations. 

The neutron measurements were performed under a 
[111] field, and IVIonte Carlo (IVIC) simulations^^^ were 
also carried out for the dipolar spin ice model^^'^^^ to 
quantify our observations. Figure 3 shows the field de- 
pendence of the magnetic intensity of the (220) Bragg 
reflection at T = Tc + 0.05 = 0.43 K. The intensity 
plateau for H < 0.8 T corresponds to the kagome ice 
state with low density monopoles. The deviation from 
the plateau as H exceeds 0.8 T indicates that monopoles 
are being created gradually, while the saturation of the 
intensity for H ^ denotes the staggered monopole 
state (Fig. 1(e)). In Fig. 3 we also show the Bragg inten- 
sity and the density of the positively charged monopoles 
obtained by the simulation. The observation shows good 
agreement with the simulation for H < H^. On the other 
hand, above there are substantially less monopoles 
than expected from the simulation, which will be dis- 
cussed below. 



We selected T = Tc + 0.05 K and = (Fig. 3 
inset) for observation of the fluctuating high- and low- 
density monopoles. At this H^T point, we measured in- 
tensity maps in the scattering plane. An intensity map of 
the kagome ice state at T = Tc + 0.05 K and = 0.5 T 
was also measured for comparison. Figure 4 compares the 
measured and simulated intensity maps. The observed 
scattering pattern of the kagome ice state (Fig. 4(a)) is 
in excellent agreement with the simulation (Fig. 4(c)), 
showing the peaked structure-^^^ at (2/3, —2/3, 0) and the 
pinch point^^^ at (4/3, —2/3, —2/3). These structures re- 
flect the kagome ice state. 

The observed (Fig. 4(b)) and simulated (Fig. 4(d)) 
intensity maps close to the critical point show a weakened 
kagome-ice scattering pattern (by the low-density state) 
and diffuse scattering around (220) (by the high-density 
state). The observation agrees fairly well with the simu- 
lation. However the diffuse scattering is less pronounced 
for the observation. We found that this discrepancy orig- 
inated from an instrumental condition of the GPTAS 
spectrometer, which has a large vertical resolution of 

= 0.25 A-^ (fuh width at half maximum, FWHM). 
We carried out the same measurement on the BT-9 spec- 
trometer. It has a smaller vertical resolution of = 0.1 
A~-^ (FWHM), which does not affect the diffuse scat- 
tering. The resulting data are shown in Fig. 4(f), which 
are in better agreement with the simulation (Fig. 4(d)) 
around (220). An interesting point suggested by this 
resolution effect is that correlations of the high-density 
monopoles are three dimensional in space, although the 
monopoles can only move in the two dimensional layers 
(Fig. 1(d)). The three dimensional correlations are con- 
sistent with the isotropic Coulomb interaction between 
monopoles. We note that the kagome-ice scattering pat- 
tern is two dimensional in nature, and thus is not af- 
fected by the vertical resolution. 

To illustrate the high- and low-density monopoles 
yielding the scattering patterns in Figs. 4(b), 4(d), and 
4(f), two typical snapshots of the monopoles of the MC 
simulation are shown in Fig. 4(e), where we depict mag- 
netic charges at centers of the tetrahedra. Lines connect- 
ing the centers of tetrahedra form a diamond lattice, 
and magnetic charges reside on its lattice points. Re- 
gions of the low-density monopoles (where black points 
dominate) produce the kagome-ice scattering pattern, 
while those of the high-density monopoles (where red 
and blue points dominate) produce the diffuse scatter- 
ing around the Bragg reflections. These critical fluctu- 
ations between high- and low-density phases reinforce 
the proposed explanation^^ of the puzzling liquid-gas 
type critical point^^^ using the similarity argument to 
phase transitions of ionic particle systems on lattices. 
Consequently they strongly suggest existence of mag- 
netic monopoles interacting via the magnetic Coulomb 
force. Further investigations of critical phenomena, 
screening of the Coulomb interaction, and effects of the 
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Fig. 4. (Color) Intensity maps measured at T = Tc + 0.05 K in the scattering plane are shown for two field values (Fig. 3 inset): the 
kagome ice state at 0.5 T (a, c) and the fiuctuating high- and low-density monopoles at He (b, d, f). (a, b) and (f) were measured on the 
GPTAS and BT-9 spectrometers, respectively, (c, d) are simulated intensities, (e) Two snapshots of the monopoles of the simulation 
corresponding to (b, d, f) are shown on the diamond lattice, in which blue, red, and black points represents +, — and magnetic 
charges, respectively. The light green circles in (b, d, f) show the high intensity regions caused by the high-density monopoles, i.e., 
ferromagnetic fiuct nations. 



anisotropic motion of the monopoles within the kagome 
lattices are of interest. 

A question, which is not pursued in the present 
study, is how monopoles unbound by the fractionaliza- 
tion move in the kagome lattice. Comparing the ob- 
served Arrhenius law with that of a study^^^ of the dif- 
fusive motion of monopoles in spin ice state, it seems 
that the interesting temperature ranges where uncon- 
fined monopoles move diffusively are roughly T > 1.5 K 
{H = 0.5 T) and T > 0.7 K {H = 0.9 T). There may be 
another interesting issue in the discrepancy between ob- 
served Bragg intensity and the classical MC simulation 
shown in Fig. 3 {H > He). We speculate that it may 
indicate the existence of quantum mechanical effects ne- 
glected in the computation. Puzzling experimental facts 
at low T were also noticed by the slow saturation of mag- 
netization^^^ and the non-zero specific heat^^^ above He 
down to very low temperatures, T < 0.1 K. For exam- 
ple, if the double spin flips shown in Fig. 1(d) can occur 
by tunneling, monopoles (or holes in the staggered 
monopole state) might move more easily than classical 
diffusion. 

Typical elementary excitations in condensed mat- 
ter, such as acoustic phonons and magnons, are Nambu- 
Goldstone modes where a continuous symmetry is spon- 
taneously broken when the ordered state is formed. This 
contrasts with the monopoles in spin ice, which are 
point defects that can be fractionalized in the frustrated 
ground states. Such excitations are unprecedented in con- 
densed matter, and will now enable conceptually new 



emergent phenomena to be explored experimentally. 
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